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CORRELATION STUDY OF LUNAR RADII OBTAINED FROM 
PHOTOGRAPHIC AND LUNAR GRAVITATIONAL FIELD DATA 

By James R. Williams and William R. Wells 
Langley Research Center 

SUMMARY 

A study has been made to determine the correlation of local lunar radii as deter- 
mined from recent earth-based photographic studies and from lunar gravitational field 
coefficients obtained from lunar satellite tracking data. The photographic study is a 
recent local radii determination performed by the Army Map Service at 734 points over 
the visible lunar surface. The lunar gravitational field coefficients used in the analysis 
are the most recent determinations from independent studies of tracking data from Lunar 
Orbiters I to IV and the Russian satellite LUNA- 10. 

The general result of the study is that all the currently determined gravity models 
are little correlated with the visible surface features of the moon except in a few isolated 
regions. 


INTRODUCTION 

Recently, new scientific data concerning the physical properties of the moon were 
obtained. The successful performance of the Lunar Orbiter series of spacecraft has 
made available preliminary estimates of the gravitational field of the moon (refs. 1 to 4 
and unpublished data from the Langley Space Mechanics Division). At the same time, 
earth-based photographic studies of the figure of the moon have been made by the Army 
Map Service (ref. 5). It is possible to study the relation between these two independent 
data types (photographic and tracking data) by a comparison of lunar radii implied by the 
two methods. The comparison is made through a computation of the correlation coeffi- 
cients of the two independently determined sets of lunar radii in various regions of the 
lunar surface. The values of the radii in each region as determined from either photo- 
graphic measurements or a lunar gravitational model are considered to be random mea- 
surements of a representative radius for the region. The regions considered represent 
isolated maria regions and isolated mountainous regions as well as combinations of the 
two. In addition, the entire visible hemisphere is considered together with a region along 
the equator between ±12° latitude over which much of the data from the Lunar Orbiter 
satellites were obtained. 



In order to calculate the lunar radii implied by the lunar gravitational field, it was 
necessary to convert the gravitational coefficients to equivalent surface coefficients as 
outlined in reference 6. This conversion requires an assumption on the internal density 
distribution of the moon. For the purposes of this analysis, the lunar density distribution 
is assumed to be uniform throughout. 

In general, the correlation coefficient between two sets of random variables is a 
measure of the linear dependence of the random variables. For the problem treated in 
this paper, the random variables are the measured and computed radii at various points 
inside a particular region on the moon. With the assumption of a uniformly dense lunar 
interior, a high value of the correlation coefficient would indicate a definite relation of 
the surface features to the gravitational anomalies, whereas a low correlation would imply 
little dependence of these two quantities. 

SYMBOLS 

Cn,m>S n>m coefficients of lunar gravitational potential spherical harmonics (n is 
degree, m is order) 

e[x] expectation value of X 

Fn,m>Hn,m coefficients of lunar surface spherical harmonics (n is degree, m is 
order) 

h(4>,\) local deviation of surface shape from sphere 

M mass of moon 

N number of measurements 

P n>m associated Legendre function (n is degree, m is order) 

r distance from center of mass of moon to an external field point 

R 0 mean radius of moon 

Rj local radius of moon determined from earth-based photographic studies 

Ro local radius of moon determined from gravitational field 
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^l,k’ R 2,k kth values of local lunar radii 


distance from center of mass of moon to an interior point 


U 


lunar gravitational potential function 


X,Y 


random variables 


a 


angle defined in sketch 1 in appendix 


selenographic longitude of point on lunar surface 




product of universal gravitational constant and mass of moon 


p.XY’^R]R 2 covariances associated with X and Y and Rj and R 2 
Pxy’^RjR 2 correlation coefficients between X and Y and Rj and R 2 

CT X’ a Y’ a Ri>°R 2 standard deviations of X and Y and Rj and R 2 
c p selenographic latitude of point on lunar surface 

A bar over a quantity denotes the mean or expectation value of that quantity. 

GENERAL CONSIDERATIONS 


A study of the relation between two recently acquired independent data types con- 
cerning physical properties of the moon has been made. In the study it was convenient to 
relate the lunar gravitational field data type to equivalent lunar radii. This set of radii 
was then compared to the set of radii derived from earth-based photographic studies 
(ref. 5) by means of a correlation study. 

Lunar Models 

The preliminary estimates of the gravitational field considered consist primarily of 
the results of analyses of the tracking data from the Lunar Orbiter series of satellites 
(refs. 1 to 3 and unpublished data from the Langley Space Mechanics Division). The 
models in references 1 and 2 were obtained through a direct analysis of the tracking data 
from the satellites. The model in reference 2 is a refinement of the model in reference 1. 
The second of these models has additional tracking data which were obtained from a wide 
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variety of orbit configurations. The model in reference 3 and the one from the unpub- 
lished data are from independent studies which were obtained from long-period analyses 
of these tracking data. An additional model supplied by the Russian results from LUNA- 10 
as given in reference 4 is also considered. 

The standard form of the lunar gravitational potential function as suggested by the 
International Astronomical Union (ref. 7) is 

U(r = £ 

where r, (p, and A are the spherical polar coordinates of the satellite relative to the 
center of mass of the moon, and C n m and S n rn are the gravity coefficients. Numeri- 
cal values for these coefficients, for the various models, are given in table I. Dashes in 
the table imply that a particular coefficient was not determined in the referenced work. 


oo n , ,n 

+ Y, Z \r) ( C n,m COS mX + Sn > m sin mX ) P n,m( sin <#>) 


n=2 m=0 


( 1 ) 


Lunar Radii 

A local lunar radius can be calculated from gravity considerations provided the 
internal density distribution and the gravitational coefficients are known. At the present 
time, the exact degree of nonuniformity of the lunar density is not known but is expected 
to be small. For this reason, the assumption of constant density was made in this analy- 
sis. With this assumption the value for the local lunar radii is obtained from the following 
equation as given in reference 6: 


R2($,A) = R 0 


1 + 


I I 

n=2 m=0 


2n + 


-)(Cn,m cos mA + S n 


m sin m V P n,m 


(sin 0) 


( 2 ) 


where cp and A are the selenographic latitude and longitude of a local point on the lunar 
surface, and C n m and S njlT1 are the gravity coefficients discussed previously. A 
brief derivation of this result is given in the appendix. 


Correlation Coefficient 


The correlation coefficient is a measure of the interdependence of two random vari- 
ables X and Y. A value of ±1 implies that the two variables are related to each other 
linearly, whereas a value of zero generally implies that no linear relationship exists 
between the variables. The usual definition of the correlation coefficient between X and 
Y is 


P XY 


^XY 

°x ct y 


(3) 
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where 


and 


MXY = E [( x - X H Y - Y)\ 




x = e[x] 

Y = E [Y] 


(4) 


( 5 ) 


A more complete discussion of the correlation coefficient and variances and covari- 
ances of X and Y can be found in reference 8. 


Correlation Coefficient for the Lunar Radii 


In this analysis, the local lunar radii Rj and R 2 as determined from photo- 
graphic and gravitational data are treated as random variables. If R^ and R 2 are 
used in place of X and Y in the equations of the previous section, the values of 

^RlR 2 ’ a Rl’ and °R2 are 

V (Rl,k - Rl)( R 2,k - R 2 ) 

^ R 1 R 2 = z. in 

k=l 


and 


2 

N 1 
= V 1 

(Rl,k - Rl) 

CTr i 

z_ 

k=l 

N - 1 

0-2 

N 

- Y ' 

( R 2,k " r 2) 

ct R2 

Z 

k=l 

N - 1 


N 


r,- V 


k=l 


- V R 2,k 

R 2= Z IT 


k=l 


"\ 




(6) 


(7) 
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where N is the number of measured values of the radius in any particular region on 
the moon, and Rj ^ and R2 > k are the kth values of a measurement of Rj and R 2 , 
respectively. Substituting equations (6) and (7) into equation (3) gives the following 
expression for the correlation coefficient between the two radii determinations (replacing 
X and Y by Rj and Rg) : 


p RlR2 


N 

I 

k=l 


l 1 * 

^ ( R l,k - R l)( R 2,k - r 2) 


N 

I 

lk=l 


Y ( R l,k - Hi) 


n 1 / 2 


N _ 
^ ( R 2,k - r 2) 

k=l 


— 1 1/2 

- \2 


( 8 ) 


Values of the correlation coefficient as given by equation (8) were computed for the 
radii determined from the five lunar gravitational models designated models L- 1 to L- 5 
in table I, and the 734 lunar radii determined by the Army Map Service in reference 5. 

ANALYSIS AND RESULTS 

For the purposes of this analysis the visible lunar surface was divided into various 
regions which represent both maria and mountainous terrain. In addition, a region between 
±12° latitude and extending across the visible side is considered as well as the entire vis- 
ible hemisphere. For each of the regions, the correlation between the lunar radii implied 
by earth-based photographic studies and that implied by the lunar gravitational coefficients 
was computed. In this analysis, surface features are considered highly correlated with 
the local gravitational field in any region for which the magnitude of the correlation coef- 
ficient is greater than 0.50. The results of the calculation for each of the gravity models 
L-l to L-5 are given in table II for each region on the lunar surface. These regions, 
together with their corresponding correlations, are shown on the lunar maps in figures 1 
to 5. 

The numbers in the first row of table II give an indication of the general correlation 
between the radii from the two data types. These numbers represent the overall corre- 
lation coefficient for each gravity model. The highest overall correlation of 0.36, based 
on the entire 734 measurements, is obtained from model L-3 which was determined from 
Lunar Orbiter I data. The next highest overall value of 0.32 is obtained from model L-2 
which was determined from a study utilizing a large amount of tracking data and a wide 
variety of orbit configurations. However, for the region between ±12° latitude and across 
the visible hemisphere (row two in table n), model L-2 has a correlation of 0.40 while 
model L-3 has a correlation of 0.22. 
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Regions 1 to 10 are located primarily in the various maria on the visible hemi- 
sphere. Taken together, these regions represent 164 individual photographic measure- 
ments of the lunar radius with some regions having as few as five and some having as 
many as 39 measurements. In regions 2 and 3, where more than half the maria region 
measurements are made, the gravity models are poorly correlated with the surface 
features. 

Regions 11 to 14 are located primarily in various mountainous ranges on the visible 
hemisphere. These regions represent a total of 183 measurements. Of the five models 
considered, model L-3 has the most consistent correlation in all four of these mountain- 
ous regions. The correlation coefficient for model L-3 varies from a low of 0.25 in 
region 14 to a high of 0.50 in region 13. Regions 11 and 12 represent more radius mea- 
surements than do regions 13 and 14 and, consequently, serve as the better indication of 
the degree of correlation in these regions. 

CONCLUDING REMARKS 

A study has been conducted to determine the correlation between the lunar radii as 
determined from earth-based photographic studies and from the lunar gravity field deter- 
mined from lunar satellite tracking data. With the exception of a few isolated regions on 
the moon, the gravity models were poorly correlated (a correlation coefficient less than 
0.50) with the surface features represented by photographic data. Of the five gravity 
models considered, the model designated L-3 has the largest correlation (0.36) for the 
entire visible hemisphere. This low value indicates, on the average, a poor correlation 
between the visible surface features and the gravitational field. The correlation coeffi- 
cients varied considerably with the gravity model in the designated maria and mountainous 
regions. However, for all models, the correlation was generally low in each type of ter- 
rain with the exception of a few isolated values. 

At the present time, it appears likely that neither of the five preliminary gravity 
models L-l to L-5 represents exactly the true lunar gravitational field. For this reason, 
it cannot be said definitely whether the correlation coefficients obtained in this analysis 
are physically meaningful. It is interesting to note, however, that the analysis on each 
model implies little correlation of the lunar surface features (lunar radii) with the lunar 
gravity field. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 21, 1968, 

129-04-01-01-23. 
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APPENDIX 


CALCULATION OF THE LUNAR RADIUS 


The surface equation for a body of uniform density can be obtained from its gravi- 
tational potential function expressed in the form 


U(r,0,X) -if 


dM 


M ( r '2 + r 2 _ 2r'r cos 


(Al) 


where the quantities r', r, and a are illustrated in the following sketch: 

P 



The denominator of the integrand in equation (Al) can be expressed in terms of 
Legendre polynomials resulting in 


U(r , 4 >,\) = £ 


1 + M 1 r ' np ” <C0S ^ 

n=2 iV1 


dM 


(A2) 


The Legendre polynomial P n (cos a) can be further expanded in terms of spherical har- 
monics as 

n 


P n (cos a) = P n (sin 0)P n (sin 0 T ) + 2 ^ 


m=l 


(n - m)'. 


(n + m)'. 


7 P n,m(si n cp) P n , m ( sin 0 T )cos m(X - A') 


(A3) 


where cp and A are the latitude and longitude, respectively, of the external point P in 
sketch 1 and cp' and A' are the latitude and longitude of the interior point Q. The ele- 
ment of mass dM can be expressed in terms of a mean density a and the latitude and 
longitude of point Q as 


dM = crr'^cos (p ' dr' d <p' dA' 

Substitution of equation (A4) into equation (A2) results in 


(A4) 
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U(r,<M) = t 


o- V 1 pir/ 2 p27r pR2(0',X') 

+ M Z pi J_ ff/2 J Q J Q r ’ +2cos ^'Pn(cos a) dr- dV d<£- 


(A5) 


U(r,0,A) = ju. 


Equation (A5) can be integrated with respect to the variable r' to obtain 

-* 77/2 p 2rr r -,n+3 

in , Q\^n 3 

n=2 

(A6) 

The equation of the surface R 2 ($,A) can be written in terms of a mean value R 0 plus 
some small perturbation about this value as follows: 


1 + M 1 ^7^^ i 0 2, '[ R 2^'’ x '>] n+3cos *’ p n<“= <*> dV «0' 


R 2(0 A) =R 0 |l + h (0,xT| 
where |h| « 1. To a first approximation 


(A7) 


|R 2 (0,X)] n+3 = R 0 n+3 [l + (n + 3)h(0,x3 

,n+3 


(A8) 


Substitution of the value given by equation (A8) for R 2 into U in equation (A6) 
results in 


Ii 


U(r,0,x) = 


1 + ®o 3 

°° , N n 

V 1 /^o\ 

rn/ 2 r 

1 + M 

L-> n + 3 V r / ^ 
n=2 

- 7t/ 2 i 


ir- 1 (^) r 7 ! c h< *'’ v)cos * ,p “ (c ° s “> dx ’ 

n=2 _7r/2 0 


d0’ 


(A9) 


It can be shown from the orthogonal properties of Legendre polynomials that the first 
infinite series in equation (A9) vanishes. This gives 


U(r,0,X) =£ 


aR n 3 V /R 0 \ n pw/2 n2ir 

1 + Z \~j J , 2 J Q h(<^)' ,A')cos 0’P n (cos a) dA' d0' 


(A10) 


aa /„ m 

n=2 ' -7f / 2 

At this point it is convenient to expand h(0,A) in terms of spherical harmonics as follows: 

00 k 


h(0,A) = J J (F k>p cos pA + H k p sin pA^P kj p(sin cp) 
k=2 p=0 


(All) 
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where the surface coefficients F^p and Hj^p depend on the function of h. Substitu- 
tion of this result into equation (A10) together with the value for P n (cos a) as given by 
equation (A3) and integrating with respect to 0 ' and X* gives 


U(r,0,A) = jr 


3 °° n / \ n , 

1 + TT" Z X fir) (^7r)( Fn ’ m COS mA + H n,m Sin mx ) p n,m(sin 0 ) 

n=2 m=0 ' 


If R 0 is determined from the condition that 

4irR 0 3 a 


= M 


then, 


U(r,0,X) = — 


oo n / x n 


DU 11 / 

1 + Z Z W (^Tl)( Fn > m C ° S mA + Hn ’ m Sin mX ) P n,m(sin 0) 


n =2 m =0 


(A12) 

(A13) 

(A14) 


A direct comparison with the form of the potential function as given by equation (A14) with 
that given in the standard form in terms of gravity coefficients, that is 


U(r,0,x) = £ 


oo n / vii 


1 + Z Z (r) ( C n,m cos mA + S n,m sin m *) p n,m(sin 0) 


n =2 m =0 

indicates that the surface coefficients are related to the gravity coefficients as 

2 n+ 1 


(A 15) 


F _ ijii -r j. 

n,m - — 3 c n,m 


ior _ 2 n -r 1 o 
n n,m ^ 3 D n,m 




(A16) 


J 


This means that in the case of a constant density assumption the shape of the surface can 
be related to the gravitational field in the following manner: 


R2(0,A) = R 0 [l + h(0,A)J = R 0 


co n 


1 + Z Z ( Fn > m COS mX + H n,m sin m ^) p n,m( sin 0 ) 


n =2 m =0 


= R r 


oo n 


1 + 


Z Z (^ 3 + COS mX + ® n J m s * n m ^) p n,m( s ^ n 0 ) 


n =2 m =0 


(A17) 


This result agrees with that given in reference 9 in which this expression is stated but 
not derived. 
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TABLE I.- COEFFICIENTS OF PRELIMINARY ESTIMATES OF LUNAR GRAVITATIONAL FIELD 


n,ra 

Model L-l (ref. 1) 

Model L-2 (ref. 2) 

Model L-3 (Unpublished data from 
Langley Space Mechanics Division) 

Model L-4 (ref. 3) 

Model L-5 (ref. 4) 


c n,m s n,m 

c n,m s n,m 

Cn,m s n,m 

Cn,m Sn,m 

Cn,m S njm 


-2.2187 X 10“ 4 
.1250 


0.2145 x 10“ 
-.3052 


-2.0596 x 10“ 
-.1661 
.2042 
-.3773 
.3012 
.1294 
.0317 
.0798 
-.1560 
.0011 
-.0082 
-.0007 
-.5505 
-.0385 
.0342 
-.0071 
-.0008 
-.0003 


10.0080 X 10-4 I 
! -.0342 ; 


.0391 
.0072 
-.0001 
I .0011 


.0445 

1 -0.0363 x KT' 

.2852 

! .1027 

.3366 


.3809 

.1115 

.0285 

.0064 

.1532 

.0537 

-.1368 



-1.9244 x 10-4 -2.06 X 10-4 

-.0897 0.0363 X 10' 4 .157 

1 .2241 .1527 .140 


.157 

0.036 X 10* 4 

.140 

-.014 

-.363 


-.568 

i -.178 

' .118 

;-.oo7 

.333 



TABLE n.- CORRELATION COEFFICIENTS FOR 
PRELIMINARY GRAVITATIONAL MODELS 


Region 

PriR 2 

for model 

- 

N 

L-l 

L-2 

L-3 

L-4 

L-5 

Visible hemisphere 

0.23 

0.32 

0.36 

0.23 

0.14 

734 

Equatorial region 

.35 

.40 

.22 

.06 

-.17 

22 

1 

.53 

-.22 

-.35 

-.07 

.52 

14 

2 

.16 

.33 

.31 

-.12 

.43 

39 

3 

-.15 

.18 

-.10 

.22 

.29 

29 

4 

-.08 

-.05 

.04 

.10 

-.10 

6 

5 

.13 

-.30 

.25 

.60 

-.52 

19 

6 

-.39 

CO 

00 

-.20 

-.24 

-.18 

22 

7 

-.57 

.52 

.62 

-.68 

.07 

5 

8 

.07 

-.47 

.16 

.25 

-.16 

15 

9 

.07 

.22 

-.10 

.05 

-.02 

10 

10 

-.12 

-.25 

.27 

-.39 

-.51 

5 

11 

.20 

.30 

.44 

.27 

.27 

67 

12 

.24 

.42 

.41 

-.10 

.54 

84 

13 

.05 

-.01 

.50 

.16 

.07 

12 

14 

-.04 

-.12 

.25 

.12 

.25 

20 



Region 

Pr 1 B 2 

N 

Visible 

0.23 

734 

Hemisphere 


Equatorial 

0.35 

22 

Region 

1 

0.53 

14 

2 

0.16 

39 

3 

-0.15 

29 

4 

-0.08 

6 

5 

0.13 

19 

6 

-0.39 

22 

7 

-0.57 

5 

S 

0.07 

15 

9 

0.07 

10 

10 

-0.12 

5 

11 

0.20 

67 

12 

0.24 

84 

13 

0.05 

12 

1 U 

-0.04 

20 


Figure 1.- Correlations for model L-l. 
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Region 

Vz 

N 

Visible 

Hemisphere 

0.36 

734 

Equatorial 

Region 

0.22 

22 

1 

-0.35 

H 

2 

0.31 

39 

3 

-0.10 

29 

4 

0.04 

6 

5 

0.25 

19 

6 

-0.20 

22 

7 

0.62 

5 

8 

0.16 

15 ■ 

9 

-0.10 

10 ■ 

10 

0.27 

5 

11 

0.44 

67 

12 

| 0.41 

i 84 

13 

0.50 

: 12 

14 

0.25 

20 


Figure 3.- Correlations for model L-3. 




Figure 4.- Correlations for model L-4. 











